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Abstract 



n - Ising spins on a random surface represented by a matrix model is studied as a 
model of the 2D gravity coupled to matter field with the central charge c > 1. The 
magnetic field is introduced to discuss the scaling exponent A , and the value of this 
magnetic field exponent is estimated by the series expansion. 



The behavior of 2d gravity coupled to a matter field has been studied by a 
matrix model in the large N limit [1]. The analytic solution by Liouville theory is 
consistent with the result of the matrix model [2]. However, these exact results are 
all restricted to the case of the central charge of the matter field c < 1, and the 
correct expression for the string susceptibility 7 st for the c > 1 remains unsolved. 

As a matrix model for 2d gravity coupled to matter field of c > 1, a multiple 
n-Ising spin model has been introduced. The series expansion [3,4] and Monte Carlo 
simulations [5,6,7] have been performed to study the universality of the phase transi- 
tion. These results suggest the existence of the intermediate region for c > 1, where 
the phase transition remains the second order with nontrivial critical exponents. 

In this letter, we extend the perturbational analysis of the previous studies [3,4] 
by adding the external magnetic field on the Ising spin, and discuss the scaling 
relation of this multiple n-Ising model coupled to 2d gravity. 

The free energy of the one-Ising model coupled to 2d gravity in the presence of 
the external magnetic field H [8] is given by 

F = -InZ 

= -In f dM+dM-exp[-tr(^Ml + \m 2 _ + a.\l. M - ge H ' M% (1) 
-ge- H Mt)\ 

where M + and M_ are Hermitian N x N matices and they represent the vertices 
on which the up and down Ising spins are placed respectively . The quantity a is a 
coupling constant related to the temperature (5 — 1/ kT as 

a = exp (-2(3) (2) 

In the large N limit, the free energy is dominated by the planar diagrams and it 
has a scaling form near the critical cosmological constant g c , 

f = <*-^<(^> < 3 > 

We restrict ourselves to the genus zero topology. The coupling constant a has a 
critical value a c at which new universal behavior appears. For the one-Ising model, 
the critical value a c is 1/4, and the string susceptibility exponent j st changes from 
-1/2 of the pure gravity case(a ^ 1/4) to -1/3 of the Ising case at a = 1/4. The 
exponent 7 st and the magnetic field exponent A are known for the unitary (p,q) 
model [9] in which the central charge c is given by c = 1 — 6(p — q) 2 /pq for c < 1; 

2 

1st = 7 (4) 

p + q — 1 
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A = (| + g)/(p + g-l) (5) 

The Ising case corresponds to p=4 and q=3 (c=l/2), and thus we have 7^ = 
— 1/3, A = 5/6. There exists a relation between these two exponents, A = (3— 7^)/4 
for c < 1. 

We extend the one-Ising model to the multiple n-Ising model, which is described 
now by 2 n matrices [3] . Each Ising spin of n-species has up and down two states, and 
there are 2 n different spin states on the vertices of the diagrams. Each Ising spin 
has the central charge 1/2 and total central charge is simply n/2, which becomes a 
central charge of the matter field c = n/2. 

For the n=2 (c=l) case, the magnetic field exponent A may be obtained from the 
general expression of the (p,q) model in the limit of p — q + 1, p — > 00 as A = 3/4. 
We will show that indeed this value is obtained by our series expansion. 

We use the series expansion about the cosmological constant g for the fixed 
topology(genus zero), from which the exponents 7 si and A are extracted by the 
ratio method. We have already obtained the string susceptibility exponent 7 st by 
the use of the series expansion up to order eight [4] . The coefficients of this series 
of the m-th derivative of the free energy behaves as 

hrn d m F/d m H = £ c k g k (6) 



The magnetic field H is applied to one spin component and (n-1) spins remain 
free of this external magnetic field. The free energy F under this magnetic field 
is easily calculated in each Feynman diagram. We develop the series expansion up 
to order eighth including the magnetic field. For example, up to order g 3 , the free 
energy of n- Ising model with a magnetic field H becomes, 

F = -2(e H + e- H )g + [lQ(e 2H + e~ 2H + 2a 2 )(l + a 2 )"' 1 

+ 2(e 2H + e~ 2H + 2a 4 )(l + a 4 ) n -V 

- {128[e 3H + e~ 3H + {e H + e~ H ){2a 2 + a 4 )](l + 2a 2 + a 4 ) n " 1 (8) 

+ ^[e^ + e -^ + ( e ^ + e ^)(3a 2 )](l + 3a 2 r 1 (8) 

+ U[e 3H + e- 3H + (e H + e~ H ){a 2 + 2a 4 )] (1 + a 2 + 2a 4 )™" 1 

+ ^[e 3H + e~ 3H + (e H + e"^)(3a 4 )](l + Sa 4 )"" 1 }^ 

+ ... 
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where the factor (1 — a 2 )~ 2 is absorbed in the redefinition of g. 

The quantity A is the inverse of the critical coupling g c and independent of 
the number m of the derivative of free energy. Therefore, A should be estimated 
consistently from the various derivatives of the free energy with respect to the mag- 
netic field, and the scaling relation should be checked with consistent values of the 
exponents j st and A. 

For n=l Ising case, we have an exact value of A as 50.625. This value is different 
from the exact solution of two matrix model A = 57. 6 [8, 10] due to our modification 
of the coupling constant g by a factor (1— a 2 ) 2 . We find that our series expansion also 
gives this critical value and by the method of the extrapolation of the successive ratio 
of Cfc for k — > oo, we obtain quite accurately the magnetic field exponent A = 5/6 
from the series of d 2 F/d 2 H. Once the correct value of A is obtained, the exponents 
•y st and A are accurately estimated in our problem, since the ratio of the coefficients 
is a monotonic function of the order of degree k. The process of the analysis may 
be shown by plotting T = —(A c — R k )k/A c against l/(k + 1) where k is the order 
of the coupling constant g and the ratio c^/ck-i is denoted by Rk- This value T 
approaches to the value of —3 + 7 st + mA for k — > oo. In fig.l, the analysis of the 
free energy and the second derivative are represented for n=l Ising model. 

For n=2 Ising case, there appears a logarithmic singularity. Due to this sin- 
gularity, the exponent 7 st becomes underestimated compared to the correct value 
j st = 0. However, the magnetic field exponent A, which is the difference between 
two exponents of the free energy and the second derivative of the free energy, is not 
influenced by this logarithmic behavior and it is correctly estimated by the previous 
ratio method explained before. We have used the values of a = 1/4 and A c = 54.0 
for n = 2 Ising model. By the same method presented in fig.l, we obtain A = 0.75 
which agrees with the correct value. 

For n=4 (c=2) Ising model, we estimate as A c = 59.2, a c = 0.23 and A = 0.79. 
We have the string susceptibility exponent j st = 0.04, which agrees with the previous 
value [3,4]. 

For n=6 (c=3) Ising model, we obtain A = 0.81 based on the value a c = 0.21 
and A c = 63.8. The exponent 7 s t is estimated by this method as 7,^ = 0.1 which 
is consistent with the previous result [3,4]. It may be interesting to note that the 
obtained value of 7 st of c = 3 appears within the error bar of the value of pure 3d 
Ising specific exponent a, which has been argued by various methods as 0.11. This 
coincidence is surprising since our model for c = 3 represents the bosonic random 
surface behavior contrary to the fermionic random surface, which appears in the low 
temperature expansion of 3d Isind model. If the sign factor due to the fermionic 
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character has no important effect and can be neglected, this coincidence may be 
understood [11]. 

For n=8 (c=4) Ising model, we take A c = 66.0, a c = 0.19. We have A = 0.825 
and 7 si = 0.2 for these values. 

The critival values of A c used here are consistent with the previous analysis [4]. 
The error bar may be ±0.05 for the estimated value of A. The accurate estimation 
of exponents for n > 8 becomes difficult within our eighth order expansion, since 
the fictitious or ghost second peak of 7 s t appears for a > a c as shown in the figures 
of ref. [4]. There is a tendency that both 7 st and A are increasing with the increasing 
value of c. It may be interesting to find the value of c, at which the branched 
polymer behavior can be seen. 

In the following, we discuss the critical behavior near a c which appears in the 
large system size, i.e. in the large k limit of the g k expansion. The coefficient c& of 
the free enegy in the large k limit becomes the partition function of n-Ising spin on 
a random surface, and we denote the logarithm of devided by k as /, 

/ = - lim (lnc fc )/fc 

~ -kiA(a,H) (9) 

This free energy / has a singulality at a c and following the conventional notations 
of the critical exponents of the specific heat a and the susceptibility 7, we have 

f~(a-a c ) 2 - a (10) 

X = d 2 f/d 2 H~(a-a c )-" (11) 
The cosmological constant A(a, H = 0) behaves near the critical value of a c as 

A(a,0) = A(a c ,0) + C(a-a c ) 1/(1, (12) 

where is a crossover exponent and it describes the shift of the cosmological constant 
due to the change of a. The crossover occurs from the pure gravity behavior to n- 
Ising behavior at a c . The second derivative of the free energy by a becomes from 
(9) and (12), 

d 2 f/d 2 a~(a-a c )^-\ (13) 
a = 2 - i (14) 
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The exponent is represented conventionally by du, where d is Hausdorff dimen- 
sion and v is a critical exponent of the correlation length. Since the specific heat 
exponent at is written conventionally by a = 2 — du, we find 



l/ct) = dv (15) 
The free energy F has a scaling form for small H and (a — a c ), 

F = to -*>^<^<F^> (16) 

Therefore, we have a scaling relation between two quantities, 

la-ad 1 /* ~ H 1 /* (17) 
and the magnetization M and the susceptibility x are expressed by 



l-A 



M = df/dH ~ ~ (a c — a)~ (18) 

x = d 2 f /d 2 H ~ H^~ 2 ~ (a-ac) 1 -^ (19) 

since we have A(a, 0) = A(a, H) + CH^ A and / is given by (9). Thus we find the 
expressions for the critical exponents a, (3 and 7 in terms of A and <fi which are valid 
for c > 1 case: 

a = 2 .^ = lzA, 1 = ^l (20) 

From these equations, it is possible to compare our estimation of the critical 
exponents with the Monte Carlo simulations, especially about (3/du and 'j/du since 
these quantities are described by the magnetic field exponent A from (15). We note 
that the previous Monte Cairo results [5,7] for the exponents (3/du, ^jdv of c = 1 
does not agree with a correct result. 

For c < 1, we have (3 = 1/2 and A = (3 — 7 s t)/4. Thus the exponent a becomes 
a = 2 r y st /{l + 7 st ) [12]. This relation does not hold for c > 1 since we have already 
found that the obtained value of A does not satisfy the relation A = (3 — 7 s t)/4. 

The magnetization M = —d\nA(a,H)/dH is obtained from the calculation of 
the inverse cosmological constant A (a, H) in our series expansion up to order eight. 
Indeed we find that the spontaneous magnetization appears below a c . However it is 
difficult to extract the accurate value of the exponent (3 or 7 from this spontaneous 
magnetization. We have obtained rather accurately the value of f3/dv — 1 — A and 
j/du = 2 A — 1 through the study of A, but the value of dv remains unknown. 
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If the spontaneous magnetization is analyzed more precisely by the higher order 
expansion, it may be possible to extract the value of = 1/ dv. 

We have also tried the specific heat analysis for n=l Ising model based upon the 
series expansion up to order ten. From the coefficient of the series expansion about 
g of the free energy F without a magnetic field, the ratio Rk = Ck/c^-i is obtained. 
Then the quantity A in (7) is estimated by a ratio method as Ak = (k+l)Rk—kRk-i- 
The second derivative of this quantity by a in (2) is proportional to the specific heat. 
For finite k, there appears a maximum peak around a critical value a c = 0.25 as 
shown in fig. 2. The extrapolation of A'j, for k — > oo is simply done by plotting this 
value against 1/k. As shown in fig. 2, we have a cusp singularity for the specific heat 
with an exponent a — — 1 and a finite size scaling relation A"^ — A" k ~ k a . This 
method seems powerful to extract the exponent a = 2 — dv, however for c > 1, 
the ghost singularity observed in ref.[4] makes the analysis to be difficult without 
extending the series expansion higher. 

As a summary, we have obtained numerically the value of the scaling exponent 
of the magnetic field A for c > 1 by the series analysis of n-Ising model, and we 
have shown how this exponent is related to other exponents. 

The author thank E. Brezin, V. Kazakov and J. Zinn- Justin for the valuable 
discussions on the present subject. He acknowledges a hospitality of Ecole Normale 
Superieure. This work is supported by the Grant-in-Aid for Scientific Research by 
the Ministry of Education, Science and Culture. 
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figure caption 

fig.l The dots show the values of T = —(A c — R k )k/A c plotted against l/(k + 1), 
where Rk are ratios of the coefficients in the series of the free energy (upper 
line) of n — 1 Ising model and of the second derivative by a magnetic field 
(lower line). The intersections of two lines with the vertical axis shows the 
exponents 3 — 7 st and 3 — 7 st — 2A. The value of A is obtained as 5/6. 

fig. 2 The second derivatives of the inverse cosmological constant Ak(a) by a = 
exp(— 2/3) are shown for k — 6, . . . , 10 in the n — 1 Ising model. The broken 
line is the extrapolation for k — > oo and it shows a cusp singularity at a = 1/4. 
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